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Context. Globular clusters are the simplest stellar systems in which structural parameters are found to correlate with their masses and lumi- 
nosities. 

Aims. In order to investigate whether the brightest globular clusters in the giant elliptical galaxies are similar to the less luminous globular 
O clusters like those found in Local Group galaxies, we study the velocity dispersion and structural parameter correlations of a sample of bright 



globular clusters in the nearest giant elliptical galaxy NGC 5128 (Centaurus A). 



■ Methods. UVES echelle spectrograph on the ESO Very Large Telescope (VLT) was used to obtain high resolution spectra of 23 bright globular 
clusters in NGC 5128, and 10 clusters were observed with EMMI in echelle mode with the ESO New Technology Telescope. The two datasets 
, have 5 clusters in common, while one cluster observed with UVES had too low signal-to-noise (S/N). Hence the total number of clusters anal- 
t ■ ysed in this work is 27, more than doubling the previously known sample. Their spectra were cross-correlated with template spectra to measure 
r> ' the central velocity dispersion for each target. The structural parameters were either taken from the existing literature, or in cases where this was 
£2 . not available, we have derived them from our VLT FORS 1 images taken under excellent seeing conditions, using the ISHAPE software. The 
velocity dispersion and structural parameter measurements were used to obtain masses and mass-to-luminosity ratios (M/L v ) of 22 clusters. 
Results. The masses of the clusters in our sample range from Mvir = 10 s - 10 7 MO and the average M/L v is 3 ± 1. The three globular clusters 
harbouring X-ray point sources are the second, third and sixth most massive in our sample. The most massive cluster, HCH99-18, is also the 
brightest and the largest in size. It has the mass (M rjr = 1.4 x 10 7 Mo) an order of magnitude larger than the most massive clusters in the Local 
Group, and a high M/L v ratio (4.7 ± 1.2). We discuss briefly possible formation scenarios for this object. 

Conclusions. The correlations of structural parameters, velocity dispersion, masses and M/L v for the bright globular clusters in NGC 5128 
extend the properties established for the most massive Local Group clusters towards those characteristic of dwarf elliptical galaxy nuclei and 
Ultra Compact Dwarfs (UCDs). The detection of the mass-radius and the mass-M/Ly relations for the globular clusters with masses greater 
than ~ 2 x 10 6 Mo provides the missing link between "normal" old globular clusters, young massive clusters, and evolved objects like UCDs. 

Key words. Galaxies: elliptical and lenticular, cD - Galaxies: Individual: NGC 5128 - Galaxies: star clusters 



1. Introduction 

The properties of globular clusters and the observed corre- 
lations between their various internal structural and dynami- 
cal parameters offer empirical constraints not only on the for- 
mation of globular clusters themselves, but also on the his- 
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tory of the host galaxy. A large number of empirical rela- 
tions between various properties, core and half-light radii, sur- 
face brightnesses, velocity dispersions, concentrations, lumi- 
nosities, metallicities, e tc., of the Milky Way globu lar clusters 
have been found (e.g. iDiorgovski & Mevlanlll994l) . Many of 
them are mutually dependent due to the fact that globular clus- 
ters have remarkably simple structures that can be reaso nably 
wel l approximated by i sotropic, single-mass lKingl (11966I) mod- 
els. iMcLauj^lm (2000) has shown that the Milky Way clusters 
are confined to the fundamental plane well defined by 2 empiri- 
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cal relations: M/L=const. and E/, ~ L 2 05 , where M is the mass, 
L the luminosity, and the binding energy of the cluster. 

The brighte st and the most mas sive globular cluster of our 
Galaxy, u> Cen (IMevlan et al.ll 1995b . is peculiar in many of its 
charact eristics: e.g. it is the most flattened Galactic globular 
cluster ( White & Shawlll987 ) and it shows strong variations in 
nearly all element abundances (e.g. iNorris & Da Costal 1 1995t 
Pancino et al.ll2002l) . A scenario that may explain some of its 
characteristics is that co Cen is the nucleus of a former dwarf 



ellipti cal galax y (Zinnecker e t alJll988tlHughes & Waller stein 
20001: iHilker & Richtlerl 12000). The same scenario was pro- 
posed for M54, another very massive globular cluster. It is 
a candidate for the nucleus of the Sagittarius dw arf (e.g. 
Bassi no & Muzziolll995l:lLavden & SaraiedinfeOOOl) . a galaxy 
that is currently being accrete d by the Milky Way. M 31 has 
4 globular clusters for which iDiorgovski et al.1 d 19971) mea- 
sured velocity dispersions cr > 20 km s _1 implying masses 
at least as large as the one of a> Cen. The most massive of 
them, Gl, shares also other particular prop erties of a) Cen, lik e 
the flattening and metallic ity dispersion dMevlan et al] 1200 II) . 
From the recent work of M a et al. (2006), the most luminous 
M31 globular cluster, 037-B327, has been suggested to be 
the most massive Local Group cluster, with a total mass of 
(3.0 + 0.5) x 10 7 M0, determined photometrically. These au- 
thors estimate the (one-dimensional) velocity dispers ion for 
037-B 327 of (72+13 km s~'). However, a later paper bv lCohen 
(2006) challenged this result, based on the measured velocity 
dispersion of cr = 21.3 ± 0.4 km s" 1 , which is comparable t o 



that of Gl (cr = 25.1 + 0.3 km s" 1 : Djorgovski et al.lll997l) . 
She concluded that 037-B327 is not the most massive cluster 
in the Local Group, and that probably M31 clusters Gl, G78 
and G280 are more massive than 037-B327. Going to galax- 
ies beyond the Local Group, ve ry similar to Gl in M31 is the 



cluster n 1023- 13 in NGC 1023 (|Larsen||2001) 



With luminosities and masses larger than globular clus- 
ters are the so-called ultra compact dwarfs (UCDs) or dwarf- 



and Virgo galaxv clusters (IHilker et al. 


1999: Drinkwateret al. 


2000: Hasesan et al. 2005t Iffilker et al. 


2007). While their ori- 



gin and relation to globular clusters is still debated in the 
literature, it has been established that very massive young 
star clusters can form in major star forming events. Such 
clusters, with masses of the order of > 10 6 , or even > 
10 7 MO dMaraston et al.ll2004trBastian et all2006l) . show sim- 
ilar scaling relations as UCDs/DGTOs, but might be differ- 
ent from the less massive globular clusters based on examina- 
tion of their mass- velocit y dispersion and mass-radius relations 
dKissler-Patig et al . 2006). -lowever, we point out that the ages 
of UCDs/DGTOs are similar to those of globular clusters, in 
contrast to the massive young star clusters forming in mergers. 

To populate the transition region between "normal" globu- 
lar clusters and more massive DGTOs, it is of interest to look 
at the massive elliptical galaxies which harbour globular cluster 
systems which are an order of magnitude more populous than 
those of the Local Group spiral galaxies. The nearest easily ob- 
servable elliptical galaxy is NGC 5 128. It has a large number of 
bright globular clusters with luminosities exceeding the bright- 
est Local Group globulars. This makes it an ideal target. 



The most recent distance determination to this galaxy is 
3.42 + 0.18 + 0.25 Mpc (the first is random and the sec- 
ond systematic error) , obtained using Cepheid PL relation 
dFerrarese et alj l2006l). H ere we use the distance of 3.84 + 
0.35 Mpc dRejkuba|[2004), which is the same value used in a 
previous work by iMartini & Hoi ( 2004 ) who presented veloc- 
ity dispersions and mass-to-light ratios for 14 bright globular 
clusters in NGC 5128. In this work we present new high reso- 
lution spectra and derive M/L ratios, thus more than doubling 
the sample of bright globular clusters with similar data in the 
literature. 

A decade ago Dubath ( 1994) presented at a conference the 
first measurements of velocity dispersions of 10 bright globular 
clusters in NGC 5128. Since these results have not been pub- 
lished in a refereed journal yet, they are included here along 
with the more recent observations of 23 clusters from UVES 
high-resolution echelle spectrograph of ESO Kueyen (UT2) 
telescope of Very Large Telescope (VLT). 

This paper is organized as follows: Sect. [2] describes the 
observations and data reduction, Sect. l3.1l shows the results of 
the cross-correlation technique for radial velocity and metallic- 
ity standard stars, while Sect. l3.2l presents the results from the 
radial velocity and core velocity dispersion measurements of 
globular clusters in NGC 5128. The comparison with previous 
measurements of clusters' radial velocity and velocity disper- 
sion is in Sect. l3.2.TI In Sect.Hthe structural parameters for 22 
clusters are presented. For those clusters which had no previ- 
ous determinations of structural parameters in the literature we 
derive them from our high resol ution ground-base d images fit- 
ting t he King profile dKingl 19621) using ISHAPE dLarsenl 19991 
200 ll) programme. In Sect. [5] we derive mass-to-luminosity ra- 
tios for the clusters and in Sect. [7] discuss the correlations and 
fundamental plane. Finally, in Sect. [8] we summarize our re- 
sults. 



2. Sample selection and observations 

Th e observations of 10 brigh t clusters (selected from the lists 
of van den Bergh et alJll98U iHesser et al. 1986; lHarris et al 



19921) wer e taken in March- A pril 1993 with the echelle mode 
of EMMI dDekker et alJll986l) . the multi-mode instrument of 
the ESO New Technology Telescope (NT T). These data have 
previously been presented at a conference dDubathlll994T) . and 
are published here together with the observ ations of 23 clusters 
obtained with UVES echelle spectrograph ( Dekk er et al . 2000) 
of the ESO VLT in April 2002. There are 5 clusters observed 
with both instruments and these have been used to check for 
the systematics in the data and errors. 

The sample of globular clusters selected for observations 
with UVES contains the brightest NGC 5128 clusters with ei- 
ther membership confirmed through published radial veloci 



ties dvan den Bergh et alll98ll:lHesser et aTll 1 986t [Harris et al 



1992L cluster names starting with VHH81, HHH86, and 
HGHH92) or the structural parameters and colours typical 



for globular cl usters in the Milky Way ( Holland et al J 1 1 999c 



Rejkuba 2001, cluster names starting with HCH99 and with 
R, respectively). 
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2.1. EMMI spectroscopy 

The first high-resolution integrated-light spectra of bright glob- 
ular clusters in NGC 5128 were obtained with EMMI at ESO 
NTT telescope during three nights, March 31 to April 2 1993. 
The red arm of EMMI was used in Echelle mode (REMD) with 
grating #10 and grism #3 (CD2), yielding the resolving power 
of 30,000, corresponding to 10 km s , and the wavelength cov- 
erage was from 4500 to 9000 A, divided among 65 useful or- 
ders. 

In total 14 spectra of 10 of the brighte st globular c lusters , 



selected from the c atalo gues of van den Be rgh et alj (1981), 
Hesser et alj d 1984ft . and Har ris et alj (1992), were secured. 



The ThAr calibration lamp spectra were taken before and af- 
ter each cluster spectrum. In addition, the following four K gi- 
ant radial velocity standard stars were observed on each of the 
three nights: NGC 2447-s28, NGC 2447-s4, HD 171391 and 
HD 176047. All the spectra were reduced with the INTER- 
TACOS software de veloped by Queloz & Weber in Geneva 
Observatory (see e.g. lOueloz et alJll995l) . 



2.2. UVES spectroscopy 

The UVES observations were carried out on the nights of 19 
and 20 April 2002 in visitor mode. The red arm of UVES spec- 
trograph was used with the standard CD#3 setting centered on 
580 nm. It is equipped with two CCDs, covering the total wave- 
length range from 4760 A to 6840 A, with a gap of 50 A cen- 
tered on 5800 A. The slit was 1" wide, giving the resolution 
of A/AA ~ 42, 000. The sky conditions were clear and the see- 
ing varied between 076 to 1'.'3, but it stayed most of the time 
around 078. 

Globular clusters observed with UVES have V-band mag- 
nitudes ranging from 17.1 to 18.8 for 22 clusters. The faintest 
observed cluster had V = 19.44. The typical exposure times 
were 1200 sec for the brighter or 1800 sec for the fainter clus- 
ters, except for the faintest 19.4 mag cluster which was exposed 
for 2700 sec. Four clusters have been observed twice and one 
cluster three times during the two night run. The multiple expo- 
sures have been averaged to increase the signal-to-noise (S/N), 
but were also reduced independently in order to provide esti- 
mates of measurement errors. The observation log for all the 
clusters is in Table Q] where we list (1) the name, (2) the obser- 
vation date, (3) the exposure times in seconds, (4) the typical 
S/N measured on the blue side of the H a line at ~ 6550 A us- 
ing the splot IRAF task. The last column lists the V magnitude 
of the clusters taken from the literature. 

Apart from the globular cluster targets, we have observed 
17 different G and K-type giant stars with a range of metallic - 
ities (-2.6 < [Fe/H] < +0.3 dex) to be used as templates for 
cross-correlation. Some stars were observed several times, thus 
yielding a total of 28 high S/N stellar spectra. The observation 
log of the template stars is in Table [2] 

After each target spectrum, globular cluster or star, we have 
obtained the ThAr lamp spectrum at the same telescope posi- 
tion. The bias and flat-field calibration data were taken at the 
end of each night. 



Table 1. Observations log for globular clusters: observations in 
1993 were done with EMMI at NTT and in 2002 with UVES at 
Kuey en VLT The nomen clature of the clusters is following that 
of the lPeng et alj (12004al) catalogue, and the magnitudes given 
in the last column are from the same catalogue where available. 
For the clusters for which there are no measurements in that 
catalogue, we take the magnitudes from the original discovery 
publications. 
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The data reduction was done both u sing the echelle pack- 
age in IRAF dWillmarth & Barnes! 19941) and the MIDAS based 
ESO-UVES pipeline dBallester et all |2000|) , where we have 
taken care to assign the wavelength calibration spectrum taken 
after each target spectrum in order to have the highest preci- 
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Table 2. Template stars observed with UVES during the 2002 run. The columns list: (1) identifier, (2) number of observed spectra, 
(3) spectral type, (4) apparent V band magnitude from the literature, (5) metallicity from the literature, (6) radial velocity from 
the literature, (7) measured radial velocity, (8) reference for catalogue value of radial velocity, and average <x CCF measured from 
cross-correlation with all the other stars for (9) lower CCD, and (10) upper CCD (see Eq.[TJfor the definition of <Tccf)- 
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References: All the radial velocities except those from D97 (Dubath et al. 1997) are compiled from 

http://www.casleo.gov.ar/catalogue/catalogue.html, which lists references to the sources: 

N04 (Nordstrom et al. 2004) : F05 (Famaev et al. 2005) : 

COR (Udry et al. 1999, see also: http://obswww.unige.ch/ udry/std/stdcor.dat). 



sion in the wavelength calibration. Due to low S/N ratio of the 
spectra, the MIDAS pipeline was not used in optimal extraction 
mode. The final spectra were normalized using the continuum 
task in IRAF and the cosmic rays were excised by hand. After 
some tests to ensure that the pipeline results were giving the 
same results as manual reductions done within IRAF, we have 
decided to later use spectra reduced within the MIDAS pipeline 
because the different echelle orders were combined in one long 
ID spectrum per CCD, offering thus the maximum number of 
lines for cross-correlation. 

3. Cross-correlation 

In order to measure radial velocities and velocity dispersions 
of all our targets we have used cross-correlation technique 
dTonrv & Davisll979h . Slightly different implementation of the 
cross-correlation technique has been adopted for the EMMI 
and UVES spectra. The comparison of the resulting velocity 
dispersion measurements for the 5 targets in common between 
the two datasets is providing a useful check on the results ob- 
tained with these two slightly different methods. 

EMMI spectra were cross-correlated with a numerical 
mask specially designed for globular clusters. Th i s has been 
described in greater detail by iDubath et al. dl99CA Il992h and 
thus we do not repeat it here. 

The globular cluster spectra observed with UVES have 
been cross-correlated with the high signal to noise spectra of 
template radial velocity stars observed during the same run, us- 
ing the IRAF task FXCOR in the RV package. All the spectra 
were Fourier-filtered prior to cross-correlation, to remove the 



residual low frequency features arising from imperfect con- 
tinuum fitting to the spectra with combined echelle orders. 
The features at frequencies higher than the intrinsic resolu- 
tion of the spectrograph were also cut. The peak of the cross- 
correlation function (CCF) traces the radial velocity, and the 
width is a function of the velocity dispersion and of the instru- 
mental width. The latter is measured by cross-correlating the 
template stars spectra with each other. Since we have observed 
a large number of radial velocity standard stars as well as gi- 
ant star templates with a range of metallicities, we could check 
that the template miss-match does not produce spurious results. 
This is described in detail in the next section. 

3.1. Template stars 

The measured radial velocities for all the stars observed dur- 
ing the 2002 run with UVES are given in Table [2] in column 7. 
They were measured by cross-correlating each Fourier-filtered 
stellar spectrum against each other. The resulting radial veloci- 
ties for each individual spectrum were averaged and we report 
here the average radial velocity and standard deviation for each 
star. In column 6 we list the radial velocity from the literature. 
These were compiled from Dubath e t"aD ( 1997 ) and the web 
database of stellar radial velocities (see table footnote). In all 
but one case the difference between our measured radial veloci- 
ties from UVES spectra and those from the literature is smaller 
than 1 km s _1 and the measurements are consistent with the cat- 
alogue values within the errors. Since the stars were observed 
with a 1.0 arcsec slit and seeing was sometimes as good as 
0.6-0.7 arcsec, part of the error in radial velocities may come 
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Fig. 1. The points with error bars are the average raw veloc- 
ity dispersion measurements (ctccf) obtained from 7 different 
template stars whose spectra were broadened to simulate the 
different input velocity dispersion (cr,„ values). The solid line is 
from Equ.Q] which is also displayed in each panel. The two di- 
agrams are for the lower and upper CCDs of UVES which have 
different cr re f values. The inserts show the difference between 
the obtained and input true velocity dispersion as a function of 
input velocity dispersion. 



from slit centering errors. The star that shows the largest differ- 
ence, HD 103295, has less certain value of the radial velocity 
and the quoted error from the literature is evidently underes- 
timated. Leaving HD 103295 out, the average difference be- 
tween our radial velocity measurements and catalogued values 
is V R (UVES) - V R (cat) = 0.07 + 0.27 km s~' indicating that 
the systematic errors due to slit centering are not significant. 
For the cross-correlation of stars with cluster spectra we adopt 
our measured radial velocity for HD 103295, and literature val- 
ues for all the other stars. 

The projected velocity dispersions (cr) for globular clusters 
are derived from the broadening of the cluster cross-correlation 
function (CCF) produced by the Doppler line broadening 
present in the integrated-light spectra due to the random spa- 
tial motion of stars. The raw measurement (ctccf) is however 



CCF 



2 , 2 



(1) 



The cr re f value is determined for both UVES CCDs by taking 
the average value of ctqcf measurements obtained by cross- 
correlating 18 selected best template stellar spectra, belonging 
to 13 different stars, with all the other stars. This same set of 
stars is used in cross-correlation of globular cluster spectra as 
well as in all the simulations (see below). Since all these stars 
are late-type giants they are not expected to exhibit line broad- 
ening due to rotation. For 3 stars, HP 6614 1, HD 107328, and 
HD 161096, Ide Medeiros & Mavorl d 1999b provide measure- 
ments of rotational velocities which are, 1.1, 1.3, and < 1.0 
km s _1 , respectively, with uncertainties which are of the same 
order of magnitude. The fact that for all the stars CCF has sim- 
ilar width (see Table |2} implies that rotation is not a concern. 
The weighted average of the stellar CCFs are 12.6 + 0.2 km s 
for the lower and 11.8 + 0.2 km s _I for the upper CCD. The 
difference in the instrumental width reflects the different reso- 
lution of the two spectral ranges. Excluding from the average 
HD 150798, the star that shows the widest cross-correlation 
peak for both spectral ranges, does not change the average 
value of (T re f. 

Figure Q] shows the validity of the Equ.Q~]for the lower and 
upper CCDs of UVES. The solid line is from Equ. Q]and the 
points represent the average projected velocity dispersion mea- 
surements obtained from 7 different template stars whose spec- 
tra were broadened by convolving each of them with Gaussians 
with known sigma (cr,„) of 5, 10, 15, 20, 25 and 30 km s _1 . 
The velocity dispersion was then measured on these broadened 
spectra by convolving them with the 18 selected best template 
stellar spectra and averaging the resulting velocity dispersions. 
The smaller inserts in each of the panels in Fig.[T]show the dif- 
ference between the average measured and input velocity dis- 
persion as a function of the input velocity dispersion, while 
the main panel displays the averages of the raw measurements 
(ctccf) at each cr,„ value. 

We have carefully selected the regions for cross- 
correlation, avoiding strong lines such as Hp, H a , sodium re- 
gion, as well as Mgb region. When we included for example 
H a line, we noticed a strong trend of cross-correlation width as 
a function of metallicity. In the final selection, this dependence 
is not present as can be seen from Fig. [2] 

In Fig. [3] we tested the dependence of the measured ve- 
locity dispersion on S/N of the input spectra. The measured 
velocity dispersions from the most broadened spectra, with 
cr m = 30 km s and with the lowest S/N, appear to be slightly 
underestimated. However, these have, as expected higher un- 
certainty and are consistent with the input values within the 
errors. 

The observations of radial velocity standards were also se- 
cured during the 1993 EMMI run and were used to check that 
the CCF has Gaussian shape and th at its width does not depen d 
on the stellar metallicity (see also iDubath et al. I ll990ill992h . 
The average sigma of the stellar CCFs, derived from 10 mea- 
surements of K-giant radial velocity standard stars observed 
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Fig. 2. Width of the cross-correlation function for stars is plot- 
ted as a function of metallicity. The measurements on the lower 
CCD are plotted with filled and on the upper CCD with open 
circles. 




S/N 

Fig. 3. Difference between the measured and simulated velocity 
dispersion as a function of S/N of the spectra. These simula- 
tions are based on randomly selected 7 template stars whose 
spectra were broadened by convolving them with Gaussians 
with widths of 10, 20, and 30 km s _I and had S/N degraded 
to simulate noisy spectra, more similar to globular cluster tar- 
gets. (See the electronic edition of the journal for the colour 
version of this figure.) 



during that same run, is 6.2 + 0.3 km s where 0.3 km s 1 is 
the standard deviation around the mean. 



3.2. Cluster radial velocities and velocity dispersions 

The initial estimate of the radial velocity for all the clusters was 
obtained by fitting the H a line. In all but one cluster spectrum 
the line was well defined and could be fitted with a Gaussian 
profile using splot task in IRAF. Then the precise radial veloc- 
ities and velocity dispersions were measured using ficcor IRAF 
task. 



Table 3. Radial velocity measurements for globular clusters 
observed with EMMI and UVES are listed in column (4). 
For clusters with multiple observations the reported value is a 
weighted average of velocities measured from individual spec- 
tra and the combi ned spectrum. The last column is the ra- 
dial velocity from Peng et al .] (E)04a) catalogue. B — V and 
V - I colours were taken from lPeng et al.1 (I2004al) when avail- 
able, otherwise from original discovery papers. They have 
been dere ddened assuming on ly foreground reddening of E(B- 
V)=0.1 1 dSchlegel et al.ll 1998b . except for clusters observed by 
Holland et alj dl999l HCH99) which have individual redden- 
ings from that work. Reddening assumed for HGHH92-C23 is 
0.31, derived from strong interstellar NaD absorption lines (see 
text for more details). 
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Radial velocity measurements for all the clusters are listed 
in Table[3] After the identifier, de-reddened (B-V)o and (V-I)o 
colours of the targets are given. The radial velocities measured 
from our sp ectra are in column 4, and the velocities from the 
catalogue of iPeng et al. (2004a) are shown for comparison in 
the last column. For clusters with both EMMI and UVES spec- 
tra the radial velocities reported are weighted averages of both 
measurements. The individual radial velocity and velocity dis- 
persion measurements for clusters with multiple observations 
are reported in Table [4] We note that the values listed as "com- 
bined" are not the averages of individual measurements, but 
rather the measurements of the radial velocity and velocity dis- 
persion from the combined UVES spectra, constructed by aver- 
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Table 4. Radial velocities and velocity dispersions measured 
on individual spectra for cluster targets with multiple observa- 
tions. The values listed as "combined" are not the averages of 
individual measurements, but rather the measurements from the 
combined UVES spectra, constructed by averaging individual 
exposures. 
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aging individual exposures. These spectra have slightly higher 
S/N and the agreement between the values obtained from the 
individual and these combined spectra indicates the absence of 
sig nificant systematic e rrors (Fig. [3). 



Dubath et alj (11997) have made detailed numerical simula- 
tions in order to understand and estimate the statistical errors 
on their radial velocity and projected velocity dispersion <r p 
measurements obtained by applying a cross-correlation tech- 
nique to integrated-light spectra. They show that statistical er- 
rors, which can be very important for integrated-light measure- 
ments of Galactic (nearby) globular clusters, because of the 



dominance of a few bright stars, are negligible in the present 
case where sampling problems are not present, thanks to the 
larger distances of our targets. An integrated light spectrum of 
an NGC 5128 globular cluster is well approximated by the con- 
volution of the spectrum of a typical globular cluster star with 
the projected veloc ity distribution. The influence of binary stars 
is negligible (e.g. Olszewski et al.l [l996. for dSph galaxies). 

The measurements of velocity dispersions for all the clus- 
ters are given in Table [5] In the first column is the identifier, 
then we list velocity dispersion measured on UVES spectra. In 
the third column the velocity dispersions measured on EMMI 
spectra are given, w hile in the last column the results from 



Mart ini & Hoi (2004) are shown for comparison. 



In Fig.|4]we plot the spectrum of the highest S/N cluster, the 
combined three 20 min exposures of HGHH92-C7, centered 
on some of the characteristic absorption lines. Overplotted are 
broadened spectra of HD 103295 made by convolving with 
Gaussians of 5 (red), 15 (blue) and 25 km s _1 (green). The 
best fitting template is the one broadened to 25 km s in 
agreement with 24. 1 + 1 .6 km s obtained by cross-correlation 
(Table|4|. The differences between the template broadened with 
25 km s _I and the cluster spectra are shown below the spectra 
in each panel. The narrow lines at 5890.0 and 5895.9 A, cannot 
be fitted by the broadened stellar templates. This is expected, 
because they are resonance lines (Na Dl and Na D2) due to in- 
terstellar ions of Nal. The equivalent width of these lines can be 
used to con strain the interstellar exti nction towards each glob- 
ular cluster dMunari & Zwittei 1997 ). The equivalent width of 
Na Dl line in this spectrum implies E(B - V) - 0.07 mag, 
somewhat lower than the average reddening of . 1 1 mag in the 
line of sight of NGC 5128 from ISchlegel et al.1 il99^maps. 
However, we note that both our measurement and the calibra- 
tion have considerable uncertainty and that the error on the red- 
dening is probably of the same size as the derived value. 

In Fig.|5]we plot cross-correlation functions for all the clus- 
ters. In each panel, next to the cluster name, the velocity disper- 
sion is given in parenthesis. The x-axis scale has been corrected 
for the relative velocity shift between the template star and the 
cluster so that the peak of the CCFs correspond to the heliocen- 
tric velocity of each cluster. All the CCFs show a single well 
defined Gaussian peak, except for R122, which has additional 
peaks at velocities 59 and at 1 14 km s . We have inspected the 
through-slit image taken at the start of the exposure as well as 
deep V-band images taken with FORS 1 under excellent see- 
ing conditions, but there is no sign of a spatially resolved blend 
at this position. The two additional peaks at the above given 
velocities are present when cross-correlating the cluster spec- 
trum with all the template stars and the height of these peaks 
is comparable, and sometimes larger than that corresponding 
to a velocity plausible for a cluster in NGC 5128. We note that 
the spectrum of this object also displays additional absorption 
features. We do not have any explanation for this, other than 
an unfortunate spatially unresolved blend with foreground MW 
stars, which is not uncommon given the low Galactic latitude 
(b= +19.4°)ofNGC5128. 

To calculate the mean and sigma of the velocity dispersion 
for each spectrum, we first calculate the straight mean and stan- 
dard deviation of the 18 independent measurements, each ob- 
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Fig. 4. The spectrum of HGHH92-C7 is plotted in gray showing the spectral regions around some of the prominent absorption 
line features: Mgb, Na D and Ha in the upper panel and Fe5270 and Fe5335 features in the lower panel. Overplotted are spectra 
of HD 103295 star broadened to simulate velocity dispersion of 5 (red line), 15 (blue) and 25 km s _1 (green). The best fitting 
template is the one widened by 25 km s _1 in agreement with 24.1 + 1.6 km s _1 obtained by cross-correlation (Table 0J. The 
differences (template-GC) are shown at the bottom of each panel demonstrating the goodness of the fit for the template broadened 
with cr = 25 km s . (See the electronic edition of the journal for the colour version of this figure.) 



tained by cross-correlating the cluster spectrum with a differ- 
ent template star. The error is the quadratic sum of the standard 
deviation and the uncertainty in the calibration of the instru- 
mental width for the CCF. This is done separately for the upper 
and lower UVES CCD. We combine the two averaged veloc- 
ity dispersion measurements through the weighted mean. Since 
the two chips have slightly different resolution and the calibra- 
tion of the cr re f is independently made, the uncertainty in the 
mean is calculated with t he average varian ce of the data using 
the following expression ( Bevingtonll 1969b : 



£[w;(*i-//) 2 ] 

i 

V (jv-i)2>, 



(2) 



where w, = 1 /cr 2 is the usual definition of weights and fi is the 
weighted mean of the x, measurements from the two chips. 

For clusters with multiple observations, the measurements 
of velocity dispersion from individual spectra are given in 



Table H while in Table \5\ we list the weighted average of the 
velocity dispersion from all the measurements on spectra taken 
with the given instrument. In par ticular the unc ertainty in the 
mean is calculated according to (Bevington 1969): 



2>; 



(3) 



3.2.1. Comparison with previous measurements 

In the upper panel of Fig. [6] we compar e our r adial velocity 
measurements with those of [Peng et aL (2004a, filled trian- 
gles). The average difference is 8 + 43 km s . Due to high 
resolution and the wide wavelength coverage the errors in ra- 
dial velocities of our spectra are significantly smaller in spite 
of their relatively low S/N. The agreement with the previously 
published measurements is excellent. 
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Fig. 5. The CCFs for all the clusters observed with UVES. The derived velocity dispersion is given in the parenthesis next to the 
name of each cluster. The x-axis scale has been corrected for the relative velocity shift between the template star and the cluster 
so that the peak of the CCFs correspond to the heliocentric velocity of the cluster. The x-axis displays one tick for every 100 
km s . All the CCFs show a single well defined Gaussian peak, except for that of R 122, which has additional peaks at velocities 
59 and at 1 14 km s _1 due to contamination by Galactic foreground stars. 



The comparison bet ween the veloc i ty disp ersions for the 
clusters in common with iMartini & Hoi ( 2004 ) work is shown 
in the lower panel of Fig. [6] The average difference be- 
twee n our velocit y dispe rsions from UVES spectra and those 
of IMartini & Hoi d2004 asterisks) is 0.14 + 1.93 km s"\ 
while there is a slightly larger difference, amounting to 2.6 + 
3.3 km s _1 , between the results obtained by these authors and 
our velocity dispersions obtained from EMMI spectra (open 
squares). 

For the five clusters in common between our UVES and 
EMMI datasets, the average difference is -0.4 + 3.0 km s _1 for 
radial velocities and 1.4 + 2.3 km s _1 for velocity dispersions. 
There is no trend with cluster brightness neither for radial ve- 
locity, nor for velocity dispersion residuals. 

4. Globular cluster structural parameters 

Holla nd et al. (1999) published the measurements of structural 



parameters for 21 globular cluster candidates in the inner part 



of Cen A, based on WFPC2 image s from the Hubble Space 
Telescope (HST). lHarris et all J2002h have used HST STIS un- 
filtered images to increase the number of clusters with mea- 
sured structural parameters to 43 in this galaxy. Most of the 
clusters in our sample have the structural parameters available 
from these two works. However, for 5 clusters observed with 
EMMI and 8 clusters observed with UVES such data do not 
exist in the literature. 

The globular cluster candidat es selected for the spectro- 
scopic observations from the list of Rejkuba ( 2001 ) have, how- 
ever, images taken with FORS1 imaging spectrograph at the 
ESO VLT under superb seeing conditions. We have thus used 
these images to derive the structural parameters for 8 clusters 
in our sample. Only one of them, R116 which is the same as 
HGHH92-C37, has previously measured parameters from HST 
imaging. The comparison between the derived parameters from 
the ground and the space data for this cluster (Table |6]l shows 
good agreement (but in the further analysis we use the more 
accurate results from the HST data for this cluster), and lends 
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Table 6. Globular cluster structural parameters, projected galactocentric distances, absolute magnitudes, masses, and mass-to- 
light measurements. The columns are: (1) ID, (2) core radius (r c ), (3) projected two-dimensional effective (half-light) radius 
(r e ), (4) three-dimensional half-mass radius (>/,), (5) concentration (c = log r,/r c ), (6) ellipticity (e), (7) projected galactocentric 
distance (R gc ), (8) absolute V magnitude, (9) virial mass (M n >) in 10 6 Mo, and (10) M/Ly. In the last column we list the source 
for the structural parameters. 
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confidence in the results from the King-profile fitting from 
these ground based images. Measurements of the profiles of 
the other 7 clusters observed with FORS1 are published here 
for the first time. 

The full descrip tion of the FORS 1 dataset used here is given 
bv lReikubal |200 1 ) and we will thus not repeat it. The most rel- 
evant parameters for the profile measurements are the seeing 
and the pixel scale. The seeing measured on the deep combined 
U-band images was (y.'52 and 0755 for the field 1 and 2, respec- 
tively. The V-band images had seeing of 0'.'54 and 0746, but 
unfortunately these bright globular clusters had saturated cores 
in V-band. The pixel scale is 072/pix. 

To measure the structural paramete rs we have used 
ISHAPE programme dLarsenlll999l l200ll) which models the 
light distribution of the cluster by convolving the assumed an- 
alytical model of the intrinsic luminosity profile of the cluster 
with the stellar PSF. The convolved model image is then sub- 
tracted from the observed image of the cluster and, via a x 2 
minimization algorithm, ISHAPE returns the best fitting model 
parameters and also produces the residual image which can be 
examined in detail. 

In particu lar for the modeling of these globular clusters the 
King profile dKingHl962 ) was assumed. In this model the core 
and tidal radii of a cluster are defined by the concentration pa- 
rameter c = log r,jr c . The implementation of the concentration 



parameter in ISHAPE is slightly different with its definition in 
the linear scale (C = r,/r c ). We shall call this ISHAPE con- 
centration parameter C in order to avoid confusion. Si nce C is 
the m ost uncertain of the fitted parameters in ISHAPE ( Larsenl 
20011) and its best fitting value depends on the initial guesses, 
we have run a series of measurements with the fixed C of 5, 15, 
30, 50, 75, 100, 150, 200, 250 and 300. In all cases the best fit 
(the lowest x 1 ) shows the smallest amount of residuals in the 
subtracted image. 

Table [6] lists the structural parameters for all our targets, 
except for the 6 clusters for which no high enough resolution 
optical images were available. In column 1 we list the cluster 
ID and its core radius (r c ) in column 2. Columns 3 and 4 report 
the projected two-dimensional half-light (effective) radius (r e ) 
and the three-dimensional half-mass radius (r/,), respectively. 
Concentration parameter (c) is listed in column 5 and ellip- 
ticity in column 6. In all cases cluster r adii are given i n par- 
sees assuming the distance of 3.84 Mpc (Reikuba 2004). This 
is the same distance used in the work of iMartini & H o (2004) 
allowing a straight- forward compar i son. W e note that using the 
shorter distance of IFerrarese et al.l d2006l) . the radii would be 
1 1 % smaller. 

Out of 6 clusters which have their structural parameters de- 
termined here for the first time, three have ellipticities larger 
than 0. 1 . We note however that significantly larger ellipticity 
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Table 5. Velocity dispersion measurements for the clusters ob- 
served with UVES are listed in column 2 and for those ob- 
served with EMMI in column 3. For clusters with several ob- 
servations the reported value is the weighted average of indi- 
vidual measurements (for each instrument separately). All the 
measurements from individual spectra are reported in Table [4] 
For compariso n in column 4 we list velocity dispersion mea- 
surements from Martini & Ho (2004). 
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has been derived for cluster HGHH 92-C37 (R116) from our 
FORS1 data with respect to work of lHarris et alJ (120021) . The 
difference in this particular case might be due to the location 
of the cluster close to the edge of FORS 1 field, where image 
distortions could have affected the measurement. The very high 
ellipticity of R 122 could be due to the fact that the image of this 
cluster is most probably blended with some foreground source 
(see also above). We have put ":" sign next to the ellipticity 
determinations from FORS 1 data that are more uncertain. 



5. Cluster masses and mass-to-light ratios 

The masses and the mass-to-light ratios for all the clusters with 
the available (or new) structural parameter measurements are 
given in Table [6] The listed masses a re virial masse s calculated 
using the virial theorem in the form (lSpitzerlll987l) : 
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Fig. 6. Comparison of the measurements from the literature and 
our results, as well as comparison between UVES and EMMI 
measurements for the clusters in common. In the upper panel 
we plot the difference between radial velocity and in the lower 
panel velocity dispersion differences as a function of magni- 
tude. In both panels the open circles are used to compare UVES 
and EMMI measurements. The filled trian g les are comparing 
radial velocities with those of iPeng et al. (|2004a[ P+04). In 
the lowe r panel our ve l ocity dispersions are compared with 
those of iMartini & Hoi d2004 MH04) for clusters observed 
with UVES (asterisks) and EMMI (open squares) separately. 
(See the electronic edition of the journal for the colour version 
of this figure.) 



where, assuming an isotropic velocity distribution, 3<x 2 is the 
mean square velocity of the stars and the cluster half-mass 
radius (r/,) is relat ed to the hal f-light (effective) radius (r e ) 
through r e ^ 3r h /4 dSpitzerll 1987h . 

1" slit at the distance of 3.84 Mpc dRejkuball2004 corre- 
sponds to 18.62 pc. Seeing better than 078 ensures that most 
of the light is in the slit. The central velocity dispersion <x has 
been estimated from the King profile fits convolved with the 
Gaussian of the width that reproduces the intensity profile with 
FWHM as meas ured along the spatial direction in the slit fo r 
each target (e.g. Djorgovski et al.lll997t lHasegan et al.ll2005l) . 
The corrections from the observed to central <x range from 4- 
10%, with the average correction of ~ 6% being valid for most 
of the clusters. 

Core radii of Cen A clusters are typically smaller than the 
resolution element, even for space based imagers, and thus the 
quoted errors for r c in Tableware probably underestimated. In 
addition, as stated before, the concentration parameter is also 
relatively uncertain, especially for the clusters that had struc- 
tural parameters determined with ISHAPE. Therefore we pre- 
fer to use the virial mass estimator as described above, rather 
than deriving the masses using the King model approximation, 
which would imply using central velocity dispersion , core ra- 
dius, and concentration parameter (lOueloz et al Jl995l) .The un- 
certainty in structural parameters imply also that the correc- 
tions of the observed <x to the central value of velocity disper- 
sion and to the infinite aperture <x are quite uncertain. In the 
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next section, when we plot <x , we apply a flat average correc- 
tion of 6% for all the clusters. 

In the calculation of the virial masses (Equ.0 cr is the mean 
value of velocity dispersion averaged over the whole cluster. In 
all our targets 1" slit samples the light to at least 3.5 r c , and 
in most cases beyond 6 r c . Comparing with t o Cen and 47 Tuc 
(Me ylan et alJI 19951 McLaughlin et al.ll2006h . the faint surface 
brightness and low velocity dispersion have negligible contri- 
bution beyond 3-5 r c , and consequently ~ 90 - 95% of the light 
of the clusters is in the slit. The corrections from the observed cr 
to infinite apertures estimated using the seeing convolved King 
profiles, are of the order of few percent. We prefer not to apply 
these rather uncertain corrections, but rather use the observed 
cr, since the corrections are smaller than the uncertainty on our 
other parameters. The negative error-bars for mass and M/L 
ratios in Table [6] include the maximum estimated aperture cor- 
rection. 

As expected from the comparison of the velocity dispersion 
measure ments, the derived m asses for the clusters in common 
with the Martini & Hoi ( 2004) sample are in good agreement 
with their virial masses for the targets in common. With more 
bright globular clusters our sample has a larger number of clus- 
ters with masses similar to, or larger than, the most massive 
Milky Way cluster a) Cen {M virial = 3 x 10 6 Mol Mevlan et al 
19951) and Gl in M31 (M virial = 7.3 x 10 6 Mo Mevlan etal 
200ll) . 



The mass-to-light ratios (M/Ly) are computed by dividing 
the derived masses with the V-band luminosities (Ly) which 
are calculated assuming the absolute V magnitude of the Sun 
to be My >Q = +4.83 mag: 

Ly = lQl-0MV-(m-M) v -A v -M v , o )] /gx 

The derived M/Ly ratios for our sample of globular clusters 
range from 0.1 to 5.9. However, we note that the cluster with 
the lowest M/Ly, R01-122, is the one that displays contribu- 
tion from (perhaps) stellar contaminants in its high resolution 
spectra (Fig. |5j. Thus its luminosity is expected to be over- 
estimated, which would then underestimate the M/Ly ratio. 
Excluding this cluster, the smallest Mj Lyis 1.1 a nd the average 
is < M /Ly >- 2.9 + 1.4, like also observed by lMartini & Hoi 
(2004). These authors note that this value is larger than the av- 
erage M/Ly of globular clusters in the Local Group galaxies 
and they explore the possible explanations for this, concluding 
that the difference is most probably real. Our analysis confirms 
their results. We discuss this further in Sect. [7] 

The errors of in the mass and M/Ly determinations re- 
ported in Table[6]include the errors from the velocity dispersion 
measurements and half-mass radii, and the errors due to aper- 
ture corrections, but do not include any systematic errors due to 
modeling or assumptions on reddening and distance. We note 
that, assuming a smaller dis t ance m odulus to NGC 5128 as de- 
termined bv lFerrarese et al. (2006), the M/Ly ratios would ac- 
tually increase by 12% and increase the difference with respect 
to Local Group globular clusters. 

The presence of a significant internal reddening in 
NGC 5128 would have the opposite effect. However, from 
the measurements of the equivalent widths of the interstel- 
lar NaD doublet we estimate that the extinction is con- 



sistent with very little or no internal reddening within the 
galaxy except for HGHH93-C23 whose spectra display mul- 
tiple and stronger NaD absorption lines (at different veloci- 
ties). Furthermore for the inner clusters, which are expected 
to suffer th e most dust reddening , the internal reddening values 
derived by Holland et alJ dl999h . are in all cases smaller than 
E(B-V) = 0. 14 mag, and mostly below 0. 1 mag. Unfortunately 
the uncertainty of the relation between NaD and equivalent 
widths and E(B - V) coupled with our noisy spectra, which 
yield high uncertainty in the NaD line equivalent width mea- 
surements, does not allow us to determine accurate reddening 
directly from the spectra. F or the inner clusters we adopt the de- 
reddened magnitudes from Holland et al.1 (119991) for the com- 
putation of the M/Ly ratios (Table©. For HGHH93-C23 we 
assume an additional E(B - V) — 0.2 mag due to dust internal 
to NGC 5128, while for all the other clusters only the fore- 
ground Milky Way e xtinction of Ay = 0.34 mag is assumed 
dSchlegel et a0l998l) . 



6. Special clusters 

6.1. Clusters with X-ray sources 

Kraft et alJ feOOlh. iMinniti et all d2004bl) . IPeng et alJ d2004ah . 



and IVoss & Gilfanovl (I2006I) have studied the Chandra X- 
ray point sources matching NGC5128 globular clusters. 
Minniti et alJ d2004bl) concluded that X-ray sources tend to be 
located in redder clusters (more metal-rich) and also prefer- 
entially reside in more luminous (massive) globular clusters. 
Three of the clusters in the present sample have been flagged 
as X-ray point sources based on the Chandra observations. 
These are C23 with luminosity L x = 1.10 x 10 38 , C21 with 
L x = 1.79 x 10 37 , and C7 with L x = 1.89 x 10 37 . Clusters C7 
and C23 are the second and the third most massive, and C21 
is sixth most mass i ve clust er in our sample. According to the 
data of IPeng et al.1 d2004al) . these three clusters are red, with 
V - 1 = 1.1 - 1.3, and luminous, with V = 17.9 to 17.2 (or 
M v = -10.3 to -11). The X-ray luminosities of C21 and C7 
are expected from typical low-mass X-ray binaries (LMXBs). 
The C23 luminosity puts this cluster on the bright tail of the dis- 
tribution of X-ray po int sources in NGC 5 1 28 globular clusters 
shown in Figure 5 of IMinniti et al.1 (I2004bl) . This can be due to 
the presence of a couple of LM XBs in this massive clust er, or 
to a single ultra-compact binary (IBildsten & Delovel2004 . The 
alternative explanation of a more massive accreting BH cannot 
be excluded, but it is really not demanded by the available data. 
In this respect, C23 appears to be similar to the globular clus- 
ter Bo375 of M31, which contains the brightest X-ray point 
source in a spectro scopically confirmed g lobular cluster, with 
L x = 2 - 6 x 10 38 dDi Stefano et alJl20oj) . 



6.2. The most massive cluster 

The brightest cluster of our sample is HCH99-18, with V = 
17.07. This is also the most massive cluster, with M vir 
1 .4 X 10 7 Mo, and the largest cluster in size (Table 6). It 
is appare ntly a metal-rich cl uster, which has normal infrared 
colours (IMinniti et al.l 1 19961) . It is located in the inner parts 
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Table 7. Properties of HCH99-1 8, to Cen and Gl. The data 



(1995, 2001) 



for the last two c lusters are from Me vlan et al 
and lHarrisI dl996l) . The t wo (V - /) colours of HCH99-18 are 
from HST photometry of Holland et al. ( 1999), and in paren- 
thesis the measurement from iPeng et alj (2004b) catalogue. 
The reddening towa rds this cluster is sep arated into foreground 
Galactic reddening ( Schleg el et alJ[l998l) . and in ternal redden- 
ing w ithin NGC 5128 along the line of sight dHolland et al.1 
19991) . 
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of NGC 5128, only 1.5 kpc away from the galactic nucleus, 
where reddening might be a problem. The internal reddening 
due to dust in NGC 5128 in the line of sight to this cluster 



is E( V - I) = 0.1 (i.e. E(B - V) = 0.06) mag (Holland et al 
19991) . This relatively low internal reddening is confirmed by 
the low total equivalent width of interstellar NaD lines. Its 
M/L v ratio is 4.7, higher than those of to Cen (M/L v = 2.4) 
and Gl (M/Ly = 3.6) globular clusters where the masses and 
Ml L were always ta ken from the same method, virial theorem 
dMevlan et alj|200lb . The (V - I) col our of HCH99-18 tak en 
from the original discovery work of Holland et alj d 19991) is 
0.99, si gnificantly bluer than the (V - I) = 1.50 in lPeng et al 
(2004b) catalogue. 

For a straight forward comparison of this most massive 
cluster (so far) in NGC 5128, with the to Cen and Gl we sum- 
marize t heir main prop e rties i n Table [7] Data for to Cen and Gl 



are from Mevlan et al. dl995ll200ll) and lHarrisI dl996l) . There 



is no available literature value for the spectroscopy of HCH99- 
18, but its redder colours indicate slightly higher metallicity if 
the same old age is assumed as for the other two clusters. This 
could explain the higher M/Ly value of this cluster. Its mass is 
factor of 2 larger than that of Gl, the more massive of the two 
Local Group clusters. 

As already mentioned in the introduction, one of the 
favoured formation scenarios for to Cen and Gl is be- 
ing remnant nucleus of stripped dwarf galaxy. This sce- 
nario has been invoked to explain peculiar properties of 
these massive clusters, such as the presence of chem- 
ical inhomogeneities, the high flattening, and the high 



central velocity dispersion, among others dZinnecker et al. 
1988; Hughes & Walle rsteinl koOdt Iffilker & Richflerl 12000; 



dispersion in HCH99-18, the three clusters share similar prop- 
erties, with the exception of low ellipticity. The ellipticity of 
HCH99-18 is actually similar to that of M 54 (e = 0.06), an- 
other Galactic globular cluster which is also considered to be 
a former nucleus of a dwarf galaxy. Thus it is possible that 
HCH99-18 also formed in a similar way. The alternative for- 
mation scenario to that of a stripped dwarf galaxy nucleus, 
could be through a merger of two or more y oung clusters 
dMinniti et al]|2004atlFellhauer & Kroupall2002l) . 

This object has the mass comparable to some of the most 
massive young massive clusters in galactic merger remnants, 
e.g. W3 and W30 in NGC 7252 and G 114 in NGC 1613 
dMaraston et all 12004 iBastian et ai]l2006l) . To the best of our 
knowledge it is the brightest and most massive old glob- 
ular cluster known to date within the distance of Cen A, 
and has similar properties to compact massive objects like 
DGTOs/UCDs observed in the Virgo and Fornax clusters, and 
therefore it definitely warrants further study. 

7. Discussion 

Figure [7] shows the relation between the absolute V mag- 
nitude and velocity dispersion. The sources of the data 
compiled from the literature and shown together with our 
Cen A clusters for comparison are given in the caption of 
the figure. The dashed line is the Faber-Jackson relation 
(IFaber & Jackson 1 19761) for bright ellipticals, while the solid 
line shows the best fit relation for Galac tic globular clusters 
from Mc Laughlin & van derMarelld2005l) . 

The bright globular clusters in NGC 5128 extend the glob- 
ular cluster luminosity-velocity dispersion relation towards 
brighter objects like DGTOs and nucleated dwarf elliptical 
(dE,N) galaxy nuclei. 

Figures [8] |9] and [10] display the relations between mass 
and velocity dispersion, mass-to-luminosity and effective (half- 
light) radius, respectively. Our bright clusters in NGC 5128 
(Cen A) are plotted together with typical old (age > 10 
Gyr) globular clusters from Local Group galaxies: Milky Way, 
LMC, Fornax dSph, M31 and M33. For comparison we also 
plot the transition object s between the globul ar clusters and 
dwarf galaxies: DGTO s ( Hasegan et al. 2005b and nuclei of 
dE,N dGeha et alj |2002). The solid lines are the best fit mass- 



sigma relation for globular clusters, the averag e dynamically 
determined M/L - 1.45 dMcLaughlinl [20001. and median 



Bekki & Freemanll2003l; iBekki & Chiball2004l) . While it is not 
possible to measure the presence or the absence of metallicity 



r e — 3.2 pc, which are for Galactic globular clu sters indepen- 
dent of mass (IMcLaughlin & van der Marelll2005l) . The dashed 
lines are the best fit relations for bright elliptical galaxies as 
discussed bv lHasegan et al. I (l2005l) . 

In all the above figures (f7l4TOb for comparison we also show 
the locations of the brightest Milky Way cluster to Cen and 
M31 cluster Gl. Since the masses of our clusters were ob- 
tained using the virial theorem, we plot virial masses where 
available: for to Cen, Gl, and M33 clusters. However, the litera- 
ture source for the other Milky Way clusters, as well as clusters 
in the LMC, SMC and Fornax dMcLaughlin & van der Marell 
20051) . presents only masses derived from King model fits, 
and the same is true for the masses of Virgo cluster DGTOs 
(IHasegan et al ] |2005l) . The masses of dE,N galaxy nuclei result 




Fig. 7. Central velocity dispersion versus absolute V band 
magnitude relation for globular clusters in NGC 5128 (Cen 
A; black filled circles, this work) is compared with old 
(age > 10 Gyr) clusters in Local Group galaxies as well as 
with transition objects between globular clust ers and dwarf 
galaxi es, UCDs/DGTOs from the work of lHasegan et alJ 
(2005) and nuclei of dwarf ellipticals in Virgo dGeha et al.1 
20021) The sources of data for the Local Group clusters are 
as follows: the Milky Way, Fornax dSph and LMC clus- 
ters are fromlMcLaughlin & van der Mare] ( 20051). M 3 1 from 
Dubath & Grillmair! dl997l), M33 from" Larsen etall d2002h . 
Gl cluster data are from Mevla n et al. d200ll) and to Cen 



from iMevlan et al.1 (11995b - For clarity error-bars are omitted 
for all but our data. They include the uncertainties in the 
correction from observed to central <x. The outlying cluster 
from our sample, R01-122 has its luminosity overestimated 
due to probable blend with foreground stars, as indicated 
by the arrow. The dash ed line is the Faber-Jackson relation 
(Fab er & Jackson! 1 19761) for bright ellipticals, while the solid 
line shows the best fit relation for Galac tic globular clusters 
from McLaug hlin & van der Marel ( 2005 ). (See the electronic 
edition of the journal for the colour version of this figure.) 



from dynamical modeling ( Geha et al. 2002|). When compar- 
ing these different systems we caution that lMevlan et al. ( 2001 ) 
pointed out that the masses from King model fits are twic e as 
large (however for counter example see lLarsen et al ] |2002l) . 

In Fig.[8]the departure from the scaling relation for globular 
clusters becomes evident. The brightest clusters in NGC 5128 
and dE,N nuclei occupy the same part of the diagram, which 
is shared also by some, but not all U CDs/DGTOs. However, as 
discussed by lHasegan et al.l d2005l) DGTOs from their sample 
have probably different formation mechanisms and thus differ- 
ent properties, with some being more similar to typical globu- 
lar clusters and others either stripped galactic nuclei or merged 
complexes of star clusters. 



Fig. 8. Central velocity dispersion is plotted as a function of 
mass. For the literature sources for the clusters from the Local 
Group galaxies and UCDs/DGTOs and dE,N nuclei from Virgo 
clusters (see caption of Fig. [7j). Solid line shows the virial the- 
orem for globular clusters and the dashed line is a dependence 
betwee n mass and <x for brig ht elliptical galaxies (Equ. 13 and 
12 from lHasegan et al . 2005). (See the electronic edition of the 
journal for the colour version of this figure.) 
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Fig. 9. Same as Fig. [8] but for mass-to-luminosity as a func- 
tion of mass. Solid line shows the average dynamically deter- 
mined MIL — 1.45 for a sample of Galactic globular clusters 
(McLaughlin 2000). ( See the electronic edition of the journal 
for the colour version of this figure.) 
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Fig. 10. Same as Fig. [HJ but for effective (half-light) radius as 
a function of mass. Solid line shows the median r e = 3.2 pc 
for Galactic globular clusters, while the dashed line shows the 
dependence between mass and effective (half-light) radius fo r 
bright elliptical galaxies (Equ. 12 from lHasegan et alj 120051) . 
(See the electronic edition of the journal for the colour version 
of this figure.) 



While the lower mass clusters do not show dependence of 
the M/L on the mass, a quite clear relation emerges for the 
clusters with masses larger than ~ 2 x 10 6 Mo(Fig. This 
is similar to the nuclei of dwarf galax ies, and UCDs/DGTOs 
( Geha et al. 2002 : Hasegan et alj|2005h and is also obeyed by 
the most massive clusters in the Milky Way (e.g. to Cen; framed 
blue triangle) and in M31 (e.g. Gl; red framed square). 

All Local Group globular clusters plotted in these scaling 
relation diagrams have ages in exce ss of 10 Gyr. The a ges of 
Cen A clusters are not well known. iPeng et al. (2004b) found 
that the metal-poor bright Cen A clusters have ages similar 
to those of Milky Way clusters, while the metal-rich appear 
younger with ages up to 5 Gyr. Since the M/Ly of a popu- 
lation of a given age increases with metallicity, it is in prin- 
ciple possible to construct a sample where the more massive 
globular clusters would be more metal-rich, thus having also 
higher M/L ratios. However we note that no relation between 
the [V - 7)o, nor (B - V)q colour and M/L is present neither in 
our data (Figure ITTb nor was noted by Martini & Ho (2004). 

Individual spectroscopic metallicities are not available for 
Cen A clusters in the literature. Therefore to explore the age- 
metallicity dependence on M/L ratio we use (B - V)q and 
(V - I)o colours and compare them to iMarastonl ( 2005 ) models 



in Fig urefTTl In the upper panels we plot the models for lKroupa 
(2001) initial mass function ( IMF), while the bottom panels 
have models calculated with ISalpeterl ( 1 1955b IMF. Looking 
only at the Cen A clusters (filled dots with error-bars), the mod- 
els with Salpeter IMF reproduce the range of M/L ratios indi- 



cating that the clusters with high M/L ratios are old and metal- 
rich. This result is very simil ar to that found for Fornax cluster 
UCDs bv lHilker et alj d2007f) . However, we point out that these 
models predict too high M/L ratios with respect to the MW old 
globular cluster M/Ls. Models with Kroupa IMF pass through 
the region occupied by MW g lobulars, but they s t ill imp ly too 
young ag es for most of them. iBruzual & Chariot! (120031) mod- 
els with IChabrierl d2003l) IMF have lower M/L ratios and thus 
fit better the r ange of M/L va l ues o f MW globular clusters (e.g. 
see Fig . 1 1 of lHasegan et al. I d2005l) and Fig. 1 1 of Hil ker et al 



(2007)). -lowever, these models do not cover the part of the 
plane where Cen A clusters lie. In principle, it is therefore pos- 
sible by choosing the IMF and different simple stellar popula- 
tion models to find a good solution for either Cen A clusters or 
MW clusters, but not both. 

The average colours of our sample of clusters in NGC 5 128 
are < (B - V) >= 0.81 + 0.07 mag and < (V - 7) >= 
0.94 ± 0. 10 mag. This is respectively 0. 14 and 0.07 mag redder 
than the average colours of the MW globulars that we use for 
comparison in Figure QT| Assuming that all the clusters have 
the same old age, the higher average M/Ly ratio for the clus- 
ters in Cen A with respect to those in MW can be explained in 
part with their higher metallicity. However, even excluding the 
reddest clusters from our sample, the average M/Ly ratio of the 
bright clusters in Cen A is still on average higher than that of 
"normal" globular clusters. As can be seen from Figure [9] their 
M/L ratios cover a range between those of "normal" globular 
clusters and those of UCDs and dE,N nuclei. 

We note also that the half-light radius is independent of 
the mass for the low mass clusters, while it increases with 
the mass for clusters more massive than ~ 2 x 10 6 Mo. This 
shows that, bright, massive clusters present a transition type 
of objects between typical globular clusters and more mas- 
sive DGTOs and dE,N nuclei. The implications of this finding 
for the young massive cl usters has been discussed in detail by 
Kissler-Patig et al. (2006). They have argued for the possibility 



of formation of such massive objects through early mergers of 
low mass stellar clusters, which might exp lain the emergence 



of mass-radius relation (for more details see lKissler-Patig et al 



(2006) and references therein). As an alternative speculation 
Kissler-Patig and collaborators suggest a possibility that "all 
star clusters form with a primordial mass-radius relation, but 
only the most massive clusters are able to retain it against the 
processes that would erase it". The result presented in FigurefT^ 
is consistent with the latter scenario, which should be further 
explored theoretically. 

Another parameter that might be linked to the formation 
scenario is ellipticity. The high ellipticity of a> Cen in our 
Galaxy and Gl in M31 have been frequently mentioned to- 
gether with other peculiarities shared by these two clusters, 
in the conte xt of the stripped galaxy nucleus formation mech 



anism (e.g. iBekki & Freeman 



gaiaxy nucieus rormation mecn- 
2 l2003b IBekki & Chibl |2004 . 



Large fraction of the massive clusters in M3 1 , LMC, and some 
in the MW show significant el lipticities (IGeisler & Ho dge 
19801: lHarrisl[l996HBarmby et alJl2000h . Ellipticities of Cen A 



clusters have be en compared to those of MW globular clu sters 
and discussed bvlHolland et all d 19991) jHarris et alj d2002h . and 
Gomez et al. (2006). They found a strikingly high fraction of 
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Fig. 11. The measured M/L ratios as a function of de-reddened (B - V) (left) and (V - 1) (right) colours for our sample of clusters 
in Cen A, and old globular clusters in the MW, are compared to expected values from SSP models of lMarastoii (2005) computed 
for two different IMFs, as shown in each panel. The ages of the theoretical predictions are shown on to the right of the curves. 
The model colours (open squares) from blue to red are for the follow ing metallicities: [Z/H]=-2.25, -1.35, -0.33, 0.0, +0.35 
and 0.67 dex. The colours of MW clusters are taken from Harrisl ( 1996 ) wet@ catalogue. (See the electronic edition of the journal 
for the colour version of this figure.) 



very elongated clusters among the luminous clusters in Cen A. 
Since our sample contains most of the clusters already exam- 
ined by these authors, it is not surprising to reach similar con- 
clusions. In addition to comparing the ellipticities of clusters as 
a function of luminosity we can test whether there is any depen- 
dence of ellipticity on the mass of the cluster - in our sample 
we find none. 



In a recent paper lFellhauer & Krou pa (2006) argue that the 
tidal heating during a close passage to the galactic center of a 
UCD may squeue the velocity distribution of its stars and there- 
fore lead to an overestimation of virial mass and M/L ratio. In 
these simulations the more compact and the more massive the 
object, the smaller is the effect on the measured velocity disper- 
sion. For the models similar to the UCDs in Virgo and Fornax 
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clusters, with core radii of the order of 25 pc and masses of 
10 7 Mo, only those that pass within 100-1000 pc from the cen- 
ter of galaxy can be significantly affected by tidal heating. 

In Table [6] we list in column 7 the projected galactocentric 
distance for the clusters in our sample. No correlation between 
the mass or M/L ratio and the galactocentric distance is ob- 
served. While without additional simulations it is not clear if 
and by how much tidal heating could be affecting the velocity 
dispersion measurements and thus mass determinations in the 
relatively (with respect to UCDs) compact globular clusters, we 
cannot exclude a possibility that some of the clusters have their 
masses overestimated due to inadequate assumption of virial 
equilibrium. However, given the range of masses, galactocen- 
tric distances and the relatively compact clusters, the explana- 
tion of the mass-M /Ly relation (Fig.|9]l is unlikely to be due 
to systematic overestimation of mass and M/Ly due to the tidal 
heating of the clusters. This relation might instead be connected 
to the formation mechanism. 

8. Conclusions 

We have presented an analysis of the radial velocities and ve- 
locity dispersions for 27 bright globular clusters in the nearby 
elliptical galaxy NGC 5128. For two targets we have confirmed 
here for the first time their membership in NGC 5128 through 
radial velocity measurements. Also, for 7 clusters we present 
the first measurements of their structural parameters from the 
King profile fitting to the high resolution ground based images. 

For 22 clusters we combine our new velocity dispersion 
measurements with the information on the structural parame- 
ters, either from the literature when available or from our own 
data, and use the virial theorem to derive the cluster masses. 
The masses range from 1.2 x 10 5 Mo, typical of Galactic glob- 
ular clusters, to 1.4 x 10 7 Mo, similar to more massive DGTOs 
and nuclei of dE,N galaxies. 

HCH99-18 is the brightest and the most massive cluster in 
our sample with M V(r = 1.4 x 10 7 Mo. With such a high mass 
and the M/L v ratio of 4.7 it is a candidate for being the rem- 
nant nucleus of a stripped dwarf galaxy. The alternative ex- 
planation could be th e merger of two or more young clusters 
( Min niti et al. 2004ah . To the best of our knowledge it is the 
brightest and most massive old globular cluster known to date 
within the distance of Cen A, and it shares similar properties 
with compact massive objects like DGTOs/UCDs observed in 
the Virgo and Fornax clusters. Therefore it definitely warrants 
further study. 

The most striking finding of our study is the emergence of 
the mass-radius and the mass-M /Ly relations for the bright 
clusters with masses larger than ~ 2 x 10 6 Mo. Figure [9] hints 
to the possible existence of two "populations" of globular clus- 
ters: (1) less massive ("normal") globular clusters, like the ones 
found typically in the Milky Way and M33, with M/Ly roughly 
independent of the mass, and (2) brighter more massive clus- 
ters, including our targets from Cen A as well as co Cen and 
Gl, with M/Ly ratios that seem to increase with increasing 
mass. Moreover, population 2 seems to link population 1 with 
more massive objects such as UCDs and dE,N nuclei. Figure 
[TOl although less clearly, suggests another difference between 



the population 1 with effective radius independent of mass, and 
population 2 with radius increasing with mass. 

This has been previously disc ussed for the young massive 
clusters dKissler-Patig et alj|2006l) in galactic m ergers and for 
DGTOs in the Virgo cluster dHasegan et al.ll2005r) . Our results 
indicate that the bright, massive, globular clusters associated 
with elliptical galaxies might present the missing link between 
"normal" old globular clusters associated with galaxies, young 
massive clusters formed in mergers and evolved massive ob- 
jects like UCDs (or DGTOs) associated with galaxy clusters. 
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